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Electronic cooling in Weyl and Dirac semimetals 
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Energy transfer from electrons to phonons is an important consideration in any Weyl or Dirac 
semimetal based application. In this work, we analytically calculate the cooling power of acoustic 
phonons, i.e. the energy relaxation rate of electrons which are interacting with acoustic phonons, for 
Weyl and Dirac semimetals in a variety of different situations. For cold Weyl or Dirac semimetals 
with the Fermi energy at the nodal points, we find the electronic temperature, Te, decays in time as 
a power law. In the heavily doped regime, Te decays linearly in time far away from equilibrium. In a 
heavily doped system with short-range disorder we predict the cooling power of acoustic phonons is 
drastically increased because of an enhanced energy transfer between electrons and phonons. When 
an external magnetic field is applied to an undoped system, the cooling power is linear in magnetic 
held strength and Te has square root decay in time, independent of magnetic held strength over a 
range of values. 

PACS numbers: 71.10.Pm, 03.67.Mn, 11.25.Hf 


I. INTRODUCTION 

Dirac [I, 2] and Weyl [3] semimetals have received an 
enormous amount of attention due to the wide range of 
exotic physical phenomena they are theoretically pre¬ 
dicted to host. For example, exotic edge states known 
as Fermi arcs [3-9] have recently been experimentally 
seen on the surface the Dirac semimetals NaaBi [10] and 
Cd 3 As 2 [11] after their theoretical prediction from first- 
principles calculations [12, 13]. There also is recent ex¬ 
perimental evidence of a Weyl semimetal phase in TaAs 
[14-16] and photonic crystals [17], after their theoretical 
predictions [18-20]. Weyl (Dirac) semimetals have lin¬ 
early dispersing excitations [which obey the Weyl (Dirac) 
equation, respectively] around the band touching points 
referred to as Weyl (Dirac) nodes. These nodes possess 
non-zero Berry curvature [21], which gives rise to non- 
trival momentum-space topology. Weyl semimetals also 
have many interesting topological properties, including 
the chiral magnetic effect [22, 23] and other phenomena 
associated with the chiral anomaly [24, 25]. The chiral 
magnetic effect is the separation of electric charge along 
the direction of an applied external magnetic field and 
occurs when band touching points have different ener¬ 
gies. There is recent experimental evidence for the exis¬ 
tence of the chiral magnetic effect in ZrTes [26] . The chi¬ 
ral anomaly causes the number of particles with a given 
chirality to not be conserved and occurs when external 
parallel electric and magnetic fields are applied. Dirac 
semimetals can be topologically protected by space group 
symmetries [2], however they are generally not as stable 
as Weyl semimetals. For an overview of Weyl and Dirac 
semimetals, see Refs. [27] and [28]. 

We are interested in exploring energy exchange be¬ 
tween electrons and phonons in Weyl and Dirac semimet¬ 
als when the electrons and lattice are at different tem¬ 
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peratures. Energy exchange with phonons is often the 
most dominate energy exchange mechanism in condensed 
matter systems [29]. As a result, energy transfer from 
electrons to phonons is a key issue with possible Weyl 
and Dirac semimetals based devices that take advantage 
of the topological properties or Berry curvature of Weyl 
and Dirac semimetals. Typically, to probe energy relax¬ 
ation, electrons are excited to high temperatures using 
an optical laser pulse. The hot electrons will then equili¬ 
brate through electron-phonon interactions and the elec¬ 
tronic temperature will approach the temperature of the 
lattice. As in normal metals and graphene, we assume 
electron-electron interactions rapidly thermalize the elec¬ 
trons among themselves during the relaxation process 
[30-32]. Transport measurements also provide a way to 
study electron-phonon interactions in Weyl semimetals, 
but as with graphene [33], resistivity due to electron- 
phonon scattering is expected to be smaller than the 
residual resistivity contribution that arises from disorder 
or electron-electron interactions [34, 35]. 

In this paper, we analytically study the energy trans¬ 
fer of electrons to acoustic phonons in Weyl and Dirac 
semimetals in a variety of situations. While we focus 
on acoustic phonons, we note that optical phonons will 
play a dominant role in cooling for electronic tempera¬ 
tures around and above the lowest optical branch. First 
principles calculations predict the optical branch to have 
a frequency of around 3.5 THz for BLi, a material that 
is expected to host a Weyl semimetal phase [36]. As¬ 
suming other Weyl/Dirac materials have a similar scale, 
our results should be applicable below some temperature 
range on the order of a few hundred Kelvin. The precise 
temperature range at which acoustic phonons dominate 
over optical phonons depends on the electron density and 
lattice and electron temperatures along with electron- 
phonon coupling strength. We are unaware of any experi¬ 
mental data on electron-optical phonon coupling strength 
for Dirac/Weyl materials, so we leave a discussion on the 
competition between cooling power of acoustic and opti¬ 
cal phonons for future work. We stress, however, there 
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should be a temperature (which is below the temperature 
of the lowest optical branch), at which acoustic phonons 
dominate cooling and our results apply. 

Alternatively, one could use our results to investigate 
the cross over between acoustic and optical phonon dom¬ 
inated cooling. By first experimentally observing the re¬ 
sults mentioned here for very low temperatures, one could 
then increase the electronic temperature until a change 
in the cooling properties is observed. When the chemical 
potential is at the nodal point, we find the temperature 
of the electrons decays as a power law in time over a few 
microseconds for Cd 3 As 2 , for example. These long-lived 
hot carriers (compared to a characteristic time scale of 
picoseconds in metals [37] when Tg is greater than the 
Debye temperature, TD)-important in calorimetry and 
bolometry [38]-exists as long as Te is less then the tem¬ 
perature of the optical branch, whereas in normal metals 
long lived hot carriers only exist for very low tempera¬ 
tures (Te Td) [30, 39]. 

In the highly doped limit, we find the temperature 
of the electrons decays linearly in time far from equi¬ 
librium and exponentially decays in time near equilib¬ 
rium. Motivated by recent electron cooling experiments 
on the two-dimensional analog (in some respects) of Dirac 
and Weyl semimetals, graphene, we also consider the ef¬ 
fect of short-range disorder in the heavily doped regime. 
For graphene, short-range disorder greatly increased the 
cooling power due to enhanced energy transfer between 
electrons and phonons [40, 41] and the relaxation rate 
can be controlled by varying disorder [42]. We show 
that such an enhancement of cooling power should be 
expected over a wide range of temperatures in Cd 3 As 2 
and other related materials. This result allows for a new 
three-dimensional material whose cooling properties can 
be controlled by disorder. Finally, in the presence of 
a moderate strength external magnetic field, the power 
loss of electrons is found to be linear in magnetic field 
strength and the temperature of electrons linearly de¬ 
cays as a square root in time with a rate independent of 
the external magnetic field. 

Our paper is organized as follows. In Sec. II, we calcu¬ 
late the power loss due to acoustic phonons for a single 
Weyl node. In Sec. Ill, we consider the effect of short- 
range disorder on electronic cooling. In Sec. IV, elec¬ 
tronic cooling due to acoustic phonons is investigated in 
the presence of an external magnetic field. Finally in 
Sec. V, we present our conclusions. Some technical re¬ 
sults are regulated to the appendices. 


II. SINGLE WEYL NODE 

We first consider the case of a single isotropic Weyl 
node. To generalize to N Weyl nodes, one can multiply 
the result for a single Weyl node by N (provided there 
is negligible scattering between nodes). To obtain the 
results for a single Dirac node, multiply the cooling power 
for a single Weyl node by two. Our approach follows the 


one taken in Ref. [29] for normal metals and Refs. [43] 
and [44] for graphene. The power loss, P, is given by 

( 1 ) 

k,a 

where /£* is the time-dependent Fermi distribution func¬ 
tion, E is the energy of the system, ek,a = ahvF\k\ is the 
quasiparticle dispersion relation for quasiparticles with 
wavevector fc, vp is the Fermi velocity, h is the reduced 
Planck’s constant, and a = ±1 labels the valence and 
conduction bands. The Fermi velocity has been exper¬ 
imentally found (via transport, optical measurements, 
angle-resolved photoemission spectroscopy) and theoret¬ 
ically predicted (via first principal calculations) to range 
from 1 X 10® m/s to 2 x 10® m/s [36, 45-52] in vari¬ 
ous Dirac and Weyl semimetal systems. In the case of 
Cd 3 As 2 , the Dirac cone is anisotropic [48]. We do not 
expect anisotropy to significantly alter our predictions. 

In this work, we take uf = 1 x 10®m/s, a value appro¬ 
priate for Cd 3 As 2 . Eq. (1) can be rewritten as a differen¬ 
tial equation for the electronic temperature, dtTf, = 
where Ce = is the electronic heat capacity. We 

note that both cooling power and electronic heat capac¬ 
ity scale with the number of Weyl nodes, thus the tem¬ 
poral evolution of Te will be independent of the number 
of Weyl nodes, under the assumption of negligible inter¬ 
node scattering. From Boltzmann’s equation, we have 

^* 4 “ = - E ( 4 “(i - - {ka o m ), ( 2 ) 

P,/3 

where 

Q 

( 3 ) 


is the transition rate between states \k,a) and [p,/3), 
M± = is the transition matrix element, = 


6{eT/±u}„), = Cfc „ —e„ 3 , is the Bose distri- 

k,p ^ k,p ’ ^ 

bution function evaluated at the temperature of the lat¬ 
tice Tl, ^ ^ ^ 

is the dispersion relation for the phonons, Cg is the speed 
of sound, 9 is the angle between k and p, p is the mass 
density of ions, V is the volume, D is the deformation 
potential constant, and 5^3 = 1 for intraband transitions 
and —1 for interband transitions. The deformation po¬ 
tential constant is just the electron-ion potential at zero 
wavevector [30]. 

In this paper, we take Cs = 2.3 x 10® m/s and p = 
7 X 10® kg/m^ which are the speed of sound and den¬ 
sity of Cd 3 As 2 [55] unless otherwise noted. The defor¬ 
mation potential has been estimated in Cd 3 As 2 to be 
in the 10 — 30 eV range from transport measurements 
[56] . Throughout this work, we take the deformation po¬ 
tential to be 20 eV. We ignore vertex corrections, which 
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give corrections that scale as CsjvF [30, 57], a small value 
in realistic systems. After some algebra (see App. A for 
details), we find the power loss, to lowest order in CsjvF 
and for arbitrary chemical potential, /r, referenced from 
the nodal point is 

P(/r,Te,TL) ft:!-— - keTL) 

r F 

/'CO 

X / dxx 

Jo 

where ks is Boltzmann’s constant. We now discuss some 
limits of Eq. (4). We note, for a general chemical poten¬ 
tial, to find the relaxation rate one must solve coupled 
differential equations (Eq. (1) and ^ = 0, where n is 
the electronic density which is assumed to be spatially 
uniform) to hnd the relaxation rate since the chemical 
potential has a temperature dependence. 

In the limit of ^ = 0, we have 

P'^ ^ F 

which gives (using Ce = , where E is the 

gamma function and rj is the Riemann zeta function). 


f{x - |3e^J.) + f{x + Pek-) 


(4) 


dt 


—frTHTe -Tl), 


D^k% r(6)7?(6) 
2pF^v%h^ r(4)r7(4)' 


( 6 ) 

We remind the reader that the temperature difference 
between lattice temperature and electron temperature is 
due to the system being hit with an optical pulse. This 
result agrees with the dimensional analysis of P, P, and 
Eq. (1) put forth in Ref. [43]. The cooling power at 
low temperatures is weak due to the high exponent of 
Te that appears in the cooling power. Physically, the 
weak cooling power of acoustic phonons in Weyl and 
Dirac semimetals is due to the small energy of acous¬ 
tic phonons, —ksT^, at a typical transition momentum 


of and the small density of states for electronic 

transitions. Far from equilibrium, i.e, in the limit that 
Tl < Te, we find Tf,(t) = —r, where tq = 

(l_|_m)3 olrJ-o 


and To is the initial temperature of the electrons. Tak¬ 
ing an initial electron temperature of 140 Kelvin, we find 
To = 36 X 10“® s. 

In the limit where Te ^ Tl, the electronic tempera¬ 
ture decays exponentially with a characteristic time scale, 
Tl = This should be compared to the low elec- 

tronic temperature (Te <C Td) cooling in metals. In this 
case, P oc Tg and we have similar slow cooling of the 
electronic temperature [30]. However, this slow cooling 
only happens in metals when Tg <C Tb. In contrast, 
slow electronic cooling in Weyl and Dirac semimetals ex¬ 
ist for a wide range temperatures (as long as Pg is much 
less than the temperature of the optical phonon branch, 
which is typically on the order of a few hundred Kelvin). 
We note that Cd 3 As 2 intrinsically has a large amount of 


charge carriers [51, 58]. Thus it is an open experimental 
question if this limit can be reached in Cd 3 As 2 . 

We now discuss cooling when the system is heavily 
doped. For Cd 3 As 2 , an experimentally well-established 
Dirac semimetal, has a Fermi energy of around 200 meV 
(in the heavily doped limit the chemical potential is the 
Fermi energy). Unfortunately, this energy scale is nearly 
the same as same the band inversion energy scale, which 
is about 250 meV [12], thus it is questionable whether 
the Dirac fermion description is applicable at finite tem¬ 
peratures for Ef = 200 meV [59]. However, recent ex¬ 
perimental progress has been made in tuning the Fermi 
level [51, 60] in Cd 3 As 2 and thus we believe our pre¬ 
dictions can be experimentally realized by lowering the 
Fermi level. To this end, we choose /i = 100 meV, a value 
well below band inversion energy scale. When ksT <C k, 
the maximum phonon momentum is 2hkF, where kF is 
the Fermi momentum. Thus, the maximum phonon en¬ 
ergy is given by hcskF- When the lattice temperature 
is below Tbg = , the Block-Griineisen tempera¬ 

ture [61], our approach breaks down. In this respect, the 
heavily doped case resembles the typical metallic case, 
where the quasielastic approximation fails below Tbg 
[43, 61]. Using 100 meV for the chemical potential gives 
a Block-Griineisen temperature of around 5 K. Given the 
recent experimental progress in tuning the Fermi level in 
Cd 3 As 2 , the Block-Griineisen temperature is also tunable 
parameter. When <C we can use the Sommerfeld 
expansion to evaluate the integral in Eq. (4). The cooling 
power is found to be 


P fti 


VD'^pPkB 


(Te-TL). 


(7) 


We also obtain (using Ce = ) 


dTe _ Te-TL _ 

dt Te ’ SkBh^vjp ^ ^ 


When Tl <C Te, the electronic temperature decays lin¬ 
early in time with a rate given by ■jp. Using the exper¬ 
imental values for Cd 3 As 2 we Hnd 7 p = 1.8 x 10^® K/s. 
Thus, the equilibration process is much faster for heavily 
doped systems compared to undoped systems. Gloser to 
equilibrium, Tg decays exponentially with a rate given by 
IpITl- 


III. SHORT-RANGE DISORDER 

We now consider the effects of short-range disorder on 
cooling for heavily doped Dirac or Weyl semimetals. Re¬ 
call in the absence of disorder, the momentum of the 
phonons are limited to 2hkF, and thus the phonons have 
small energies, -^k- With disorder, phonon momentum 
is no longer restricted and may reach up to kBT/cg [42]. 
This provides a boost to cooling power since the phonons 
can take away more energy from the electrons. For low 
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impurity concentrations, this process can be described 
by dressing the electron-phonon vertex. Following the 
formalism developed in Ref. [42], we derive the transi¬ 
tion matrix elements and analytically find the power loss 
(see App. B for a derivation of the transition matrix ele¬ 
ments) for disorder described by the following zero-range 
potential V{r) = d{r — rj){l + (7^)12, where rj is 
the location of the jth impurity and cr^ is the third Pauli 
matrix. This formalism has had success in describing 
electronic cooling in graphene. 

As first mentioned in Ref. [42], allowing for this type 
of generalized disorder that depends on the spinor struc¬ 
ture of the wave-function allows one to obtain a ma¬ 
trix element linear in u. In other words, pure scalar 
disorder gives a quadratic dependence on the disorder 
strength and thus the power will scale as u'^, which is 
small for weak disorder and won’t provide efficient cool¬ 
ing. Physically, this type of disorder might arise from 
differences of sublattice potential [42] or magnetic im¬ 
purities [62, 63]. The transition matrix elements are 

X which we 

plug into Eq. (1) and then perform an angular average 
to find 


^ V-Kv{^)D^k% ,rp4 _ 

p{hCs)^VFl 30 ® ^ 


(9) 


where I is the mean free path and v{p) is the density 
of states at the Fermi energy. The ratio of power loss 
for disorder to the normal momentum conserving process 
[Eq. (7)], after linearizing Eq. (9) in Tg — Tl is 

^ ^ ( Te 

P 120 kpl \Tbg 

The mean-free path in Cd 3 As 2 is on the order of a 100 nm 
[64]. We thus take kpl = 40. Disorder assisted cooling 
then dominates if Tg > 2Tbg- This result is insensitive 
to the precise value oi kpi because of the cubic root in 
Eq. (10). The ratio for arbitrary values of electron and 
lattice temperature is 



Pd Tj + T^Tl + T,Tl + Tl 

P UQkpl TIq ' ^ ’ 

Taking Tg = 50 K, = 10 K and Tbg = 5 K, we find 
that cooling power is enhanced by a factor of 250 in the 
presence of short-range disorder. The time evolution for 
Tg > Tp is 


Te{t) 


To 

"x/l “h 2(7iTo 


( 12 ) 


77^k 

where a = gop(fte )3vF{kFi) ' cooling time can be con¬ 
trolled by tuning the amount of disorder. This possibility 
was first suggested in Ref. [42] for graphene. More specif¬ 
ically, the ratio of time it takes to cool to some temper¬ 
ature (which is greater than the lattice temperature) for 


two different disorder strengths is the ratio of the mean 
free paths. 

We note that scattering from Coulomb impurities will 
likely dominate electronic transport due to inefficient 
screening of three dimensional semimetals [34, 65], how¬ 
ever we expect electronic cooling by acoustic phonons 
to be dominated by short-range disorder as in graphene 
[ 66 ]. Physically, this is due the fact that Coulomb dis¬ 
order doesn’t have any spinor structure and vanishes to 
first order in disorder strength. 


IV. ELECTRON COOLING IN A MAGNETIC 
FIELD 

In this section we consider the effect of an applied ex¬ 
ternal magnetic held on electron cooling in the cold, neu¬ 
tral limit. Most of the interesting physics of topological 
semimetals involve the presence of an external magnetic 
held [24, 25, 67]. We consider the following low-energy 
Hamiltonian of a Weyl semimetal with two nodes in an 
external magnetic held (with h = 1) [ 68 ], 



where A is the vector potential, ip and '0 = are 

four component spinors, e is the electric charge, 7 are 
the four-dimensional Dirac matrices in the chiral repre¬ 
sentation and 7 ^ = i 7 ° 7 ^ 7 ^ 7 ^. The Weyl nodes are sep¬ 
arated in momentum space by the vector 2b. We take 
the two Weyl nodes to be at the same energy, i.e. the 
zero-component of the four vector b is zero. Our results 
also apply for Dirac semimetals when 5 = 0. Recall Weyl 
nodes at different energies give rise to the chiral mag¬ 
netic effect, which generates an electrical current parallel 
to the external magnetic held. This electrical current will 
cause phonon drag [69]. Thus, the phonons will dissipate 
energy along with transporting it. While this situation is 
very interesting, it is beyond the scope of this work. The 
Weyl nodes are taken to be separated in the z-direction 
and the magnetic held, B, is also taken to lie in the z- 
direction. We leave any possible dependence of the cool¬ 
ing power on the angle between b and B for future work. 
We ignore the effects of the magnetic held on the ions 
due to their large mass. 

The cooling power can be written in terms of the imag¬ 
inary part of the self-energy as originally derived by Ko¬ 
gan [70]. The power loss is given by 

Q 

(14) 

where Imn^^((f,/3e, w) is the imaginary part of the 
phonon self-energy and Im D{q,uj) = T:(5(u} — u}q) — 5(u} + 
ojq)) is the imaginary part of the phonon Green’s function 
[71]. The phonon self-energy to one loop order is given 
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by [30], 

Tf^{q,iujn) = 

X! X! Ti-[ 7 °G(fc, ipmh°G{k + q, ipm + iuJn)], 

iPm 

(15) 

where G{k, ipm) is the Greens function of the electrons in 
the presence of an external magnetic field, 
is the electron-phonon coupling strength and /3e is the 
inverse electron temperature. We derive the imaginary 
part of the phonon self-energy in the presence of a mag¬ 
netic field in App. C. We find the exact expression for 
the imaginary part of the lowest Landau level contribu¬ 
tion to the self-energy (after analytically continuing to 
real frequencies) is, 

ImP‘'[n(g,fl)] = ^ - XvFq.)e-'^, 

(16) 


where = y ^ is the magnetic length and q‘j_ = q^+qy- 
We observe that the imaginary part of the self energy is 
independent of temperature and chemical potential. A 
similar feature is seen in the current-current correlation 
function for Weyl semimetals in the presence of an exter¬ 
nal magnetic field [68]. 

After plugging Eq. (16) into Eq. (C26), the contribu¬ 
tion of the lowest Landau level at /i = 0 at a lattice tem¬ 
perature of zero to the power loss is (restoring factors of 
h) 


Pb 


VD‘^u}% 

320TT'^VplgpC‘^ ’ 


(17) 


where wo is the Debye frequency (See App. C for some 
keys steps in this derivations). In this case, ujd acts as a 
high-energy cut off. This result is valid to lowest order in 

^ and in the limit where as = 

The limit that physically reasonable. 

For example, taking Tg = 30 Kelvin, hvp^ = 1200 
Kelvin (this magnetic energy scale can be reached with 
only 9 Telsa for vp = 10® m/s), ^ = 100, and Tp = 140 
Kelvin, we find x i —^ « .15. For Weyl and Dirac 

^ Cg 1 e 

semimetals, the Debye temperature can range from 140 
Kelvin in Au 2 Pb [35], 200 Kelvin in Cd 3 As 2 [55] to 420 
Kelvin in the pyrochlore iridates [72, 73]. While we have 
derived this result for zero lattice temperature, it is ap¬ 
plicable when Te ^ Tp. We note this result is only valid 
when the scattering between Dirac/Weyl nodes of differ¬ 
ent chirality is weak. 

The contribution from the lowest Landau level will 
dominate as long as as ^ 1- All higher Landau lev¬ 
els are exponentially suppressed by where n is 

the nth Landau level. (We note at finite chemical po¬ 
tential, apn must be greater than /x/3e to suppress the 


nth Landau level.) From this, we find (using the heat 
capacity, Gg = [74], which is valid when ^ 1) 


5Te _ 1 _ 3D^k%T^ 

dt Tf,' pvF{k,Cs)^' 


(18) 


Notably, this rate is independent of the magnetic field 
and the distance between Weyl nodes (provided scatter¬ 
ing between nodes is negligible). We note there will be 
small corrections due to the small contribution of higher 
Landau levels. The electronic temperature decays as 


Tg(t) = (19) 

Taking ujd to be 140 Kelvin and Cs = 5 x 10® m/s, we find 
7 b = 40x 10®^ K/s. We do not use the speed of sound and 
Debye temperature of Cd 3 As 2 for this calculation since 
that material has a larger ratio of Fermi velocity to speed 
of sound and Debye temperature then other Dirac/Weyl 
materials (this makes the limit ^ harder to 

reach). Furthermore, to see quantum limit transport in 
Cd 3 As 2 one needs fields of 43 Tesla [75]. This is due 
to the large Fermi surface of Cd 3 As 2 . As a result of 
such high fields and low temperatures, one might expect 
electron-phonon coupling to be modified [76]. 

Finally, we remark that these results are rather unique 
to Weyl or Dirac semimetals and one does not generically 
expect to see electronic cooling dominated by the lowest 
Landau level in normal metals. This is because the mag¬ 
netic energy scale for Weyl or Dirac semimetals is much 
larger than that of normal metals [74]. More explicitly, 
the magnetic energy scale for Weyl or Dirac semiinetal 
with vp = 10® m/s is 1200 Kelvin for a 9 Telsa magnetic 
field. For the same applied field in a normal metal it is 
12.6 X — Kelvin, where m is the effective mass and rup 
is the electron mass. For most metals, •nif.lm is on the 
order of unity [77]. 


V. CONCLUSION 

In this work, we have analytically studied the cool¬ 
ing power of acoustic phonons as a function of doping 
level, disorder, and externally applied magnetic fields. 
Our main results are in Eqs. (5), (7), (9) and (16), along 
with the corresponding decays for the electronic temper¬ 
ature, Tp, in each case. Importantly, we find disorder can 
effectively be used to control the cooling power in Cd 3 As 2 
and other closely related materials. We stress that we 
have ignored electronic cooling from optical phonons and 
that our results are only valid for some temperature be¬ 
low the temperature of the lowest optical branch. This 
temperature depends crucially on the chemical potential, 
electronic-optical phonon coupling strength, and electron 
and lattice temperatures. We leave these material spe¬ 
cific details as an open question. In future work it would 
be interesting to study the effect of Fermi arcs and Kondo 
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impurities [78] on electronic cooling, as well as interac¬ 
tions [79-83]. 

Note Added - Just prior to completion of this work, we 
noticed experimental results on the cooling by phonons 
in Cd 3 As 2 for temperatures far above the temperature 
of the lowest optical branch [84]. It was suggested in 
Ref. [84], that hot carriers and optical phonons equili¬ 
brate rapidly (500 x 10“^^ s) followed by slower cooling 
(10“^^ s) through the emission of acoustic phonons by 


the decay of optical phonons or hot carriers. 
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Appendix A: Power Loss of Single Weyl Node 


In this section, we provide some key steps in the derivation of Eq. (4), starting from Eq. (2). It is first helpful to 
seperate the power into two terms. Find and Pspom depending if they describe induced transitions or spontaneous 
transitions [44]. These two terms are given by (with h= \) 


p^nd{g,TM = 


9 pP ka 


k,p 


and 


Psponig.n) = 


9 pP ka 


k,p 


(Al) 


(A2) 


We note that Pspon{g,Te,Te) = Pind{g,Te). As such, we only need to evaluate Find [44]. We now consider the limit 
Cs <C vf as discussed in the main text. In this limit, we can neglect inter-band transitions, i.e. a ^ (3 [43]. It is 
instructive to consider each term in the a sum separately. For a = 1, using the delta function and the identity, 
/o°° deJ(e - e^+i)ff(e) = we have 

_ poo 

P“f(M,r„T,.) = -27rEE / de{e-el)wl^[f{e)-f{el)]N^{w,)S{e-eyoj,)S{e-el+^). (A3) 

9 P 

We can rewrite — vfP — Wq) as to evaulate the p integral. This gives 

PS'ip t..Tl) = /ii f de/__dx(e-( 


1 -b 


(e — cq) -I- vpqx 


ipvp (27r)® J "‘J^q J_i " 1 y \/{e- cq)^ + Vpq'^ + 2(e - cq)vFqx 

X [/(e) - /(e - cq)]N^{cq)5 - \J{e - cqY + v],q^ + 2(e - cq)vFqx'^ . 


(A4) 


After making the q integral dimensionless, we have (to lowest order in 


^OO poo 




47r(27r) 

(27r)e 


dy I de I dxe^y^ ( 1 + 


-I- y 2 -I- 2eyx 


V'e^ + + 2eya;^ . 

’ + 2eyxJ de \ J 


(A5) 

The remaining delta function can be rewritten as J e"^ + y'^ + 2eyx^ ~ A 2 e)y • Using the delta function to 

evalute the x integral, we find 


p2€ poo 


Pt=dHpi,Te,TL) = -2 Vf 


47r(27r) 
dpvf,PL (27r)6 /q Jg 


dy / dee^y^ 2 


=_E 


2e^J de pVpI3l{2tt)^ 6 Jq de' 


dee®^. (A6) 


After integrating by parts and making the e integral dimenionless, we have 


p-=\fj„T,,TL) = ^^(kET^fiksTL) / dxx^fix-P,g). 

Jo 


(A7) 
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After performing a similar calculating for ^{fj,,Te,TL), we find Find is given by 

T/r )2 1“°° f \ 

Pind{fJ-,Te,TL) = - g—^{kBTe)^{kBTL) / dxX^ [ f {x - I3efl) + f {x + j3efJ.) ) ■ 

Jo \ J 

The total power loss is then 


P=- 


VD^ 

pVpTT^ 


{kBTe)^{kBTe - PbTl) J dxX^ (^f {x - (Jep) + f {x + ■ 


(A8) 


(A9) 


Appendix B: Derivation of Transition Matrix Element 


In this section we derive the transition matrix element in the case of weak zero-range disorder described by the 
following potential 

= + (Bl) 

This derivation generalizes the one in Ref. [42] for graphene to three dimensional topological semimetals. We assume 
that the concentration of disorder is low and can dress the electron-phonon vertex by scattering off a single impurity. 
The exact transition matrix element is given by 

M± = {k’\MlG{p)f -f fG{p)Ml + fGMlG{p)f\k) (B2) 


where G{p) is the free electron Green’s function, T is the scattering operator (or T-matrix) for a single impurity. The 
scattering operator to lowest order in disorder strength is taken to be the Fourier transformed impurity potential. 
We now make some approximations of the free electron Green’s functions, similar to the ones made in Ref. [42] for 
graphene. As mentioned in the main text, this formalism has been successful in providing understanding experimental 
results of electron cooling in graphene. We expect phonons with momentum fc^T/cs, to dominate cooling. As such, 
we expect the virtual electrons to have much larger momentum than incoming and outing going electrons (fc, k' p). 
This allows one to approximate the electron Green’s function, when the virtual states have an energy hvpp ^ kBT, p, 
as G{p) = Plugging this into Eq. (B2), we find 

\M±\^ = \M\^ = X (B3) 

where no is the impurity concentration. Here, the summing over impurities is done after squaring M± [66]. For 
simplicity, we use the angular average of j(fc'j((f x which is q^/2. li p ^ kBT, we can approximate the sum 

over k and k' as f J dede'. The power loss is then 


P = n(/x)V^(^lA/+l2w,y' de/(e)(l-/(e + cc,))Arf + lM_p(-c.,) J de/(e)(l -/(e - u;,))(Arf + 1)^ . 

Evaluating the remaining e integral, we have (defining the mean free path, I = 2Tru^no’^G) 

^ VTTv{p)D^kl TT^ 4 _ 4 

pihcsfvFl 30 ^ " 


(B4) 


(B5) 


Appendix C: Phonon Self-Energy in a Magnetic Field 


In this section we derive the imaginary part of the phonon self-energy in a magnetic field. The Green’s function 
(for a given chirality, x) for ^ Weyl semimetal described the Hamiltonian in the main text is given by [68] 

(-ir 


GHu;,k,kF) = E - \vF{k. - x6)7'){iP_L„(2fcl/|) - iP+L„_i(2fcl/|)}+ 

1 


A=± n=0 


2\vF{kF ■ 7J-)ii-i(2fciZ|)) 


U) -\- IJ, — XEn 


(Cl) 










where L“ are the generalized Laguerre polynomials, 'P± 


i ^1 ± isign(ei?) 7 ^ 7 ^^ and 


En = VF\ {kz - xW + 2n 


\eB\ 


(C2) 


Following Ref. [85], we rewrite our Green’s function in a mix of real-space and momentum space coordinates. The 
partial Fourier transform of the Green’s function is given by 




(27r)2 


The inverse partial Fourier transform is 


G{uj,kz,k±) = ^ J drie fc^,rj_). 


(G3) 


(G4) 


The hybrid real-space/momentum-space Green’s function is then 


GXK£.rU = .|:2je-4 


EE 

A=± n—0 


(- 1 )" 


(E^lo - AuF(fcz - xbh^){E-Lr, 


2/2 


— V^Ln-l 


0/2 

ALr 


}- 




-\El' 


(G5) 


The total Green’s function for both chiralities is then 

G(cG,fc„r-l) = ^ G^{u:,kz,r^)V^. (G6) 

It is convenient to introduce the spectral function 

A{Lo,kz,r:0 = ^(G^=o(w-ie, fc2,r_L)-G^=o(w-l-ie, fcz,rj_)) = ^ kz,r±)V^, (C7) 

x=± 


where 


.U 1 -- 


A=± n—0 


EE 


i-iy 


{E^lo - ^vpikz - xb)l^){E-Lr, 


2/2 

ZIb 


— E+Ln-1 


2/2 

ZLb 


}- 


2z^Xir± ■ - Ai?^), (G8) 


as done in Ref. [68]. The spectral function and Green’s function are related by 

A /'°° da;.4(u;, fc^, rj_) 

G{iuJn,kz,r±) = / —-^-. 

7-00 lUJn+ H-UJ 

As discussed in the main text, the phonon self-energy to one-loop order is given by 

q, iu)n) = ^ ^ Tr[ 7 °G(fc, ipm)j°G{k + q, ipm + iw„)], 

'^'Prn ^ 


(G9) 


(CIO) 


where the trace is over spinor indicies. After performing the Matsubara sum and analytically continuing (iun —> fl-|-i? 7 ) 

k 


n {B,q,n) = —— / dw / dw- - —^- > Tr 7 ’’A(fc,a;) 7 M /c-tg,w . 

V J J UJ — UJ' — \l — IT] ^ 


(Cll) 
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We are only concerned with the imaginary part of the phonon self-energy. Using the identity = P(^) -I- iTT5{a) 
(when a is real), we have 


Imll^ {B, q,Vt) = g / dw ne(w — g) — ne{uj — — g) 


dk, f d^/cj 


Tr[y° A{k,ui)j^A{k + q,uj — id)]. (C12) 


27r J (27r)2 

Switching to the real-space/momentum space spectral function via Eq. C4, we have 

lmlf’^{B,q,n) = J du;(^ne{u}-g)-ne{uj-n-g)'^ j j . 


There are two different real-space integrals that need to be evaluated. They are as follows 




and 


=27r/V(-l)("+"')e-^Ll-’ 


„2 ,2 / /2 

- -K, ^ T n' —n ( ^A.^B \ r n—n' 


Ly 


qVl 


(CIS) 

(C14) 

(C15) 


After performing the trace, we are left with three terms that group by their Laguerre polynomials and the imaginary 
phonon self-energy can be written as sum of two terms, Ii + l 2 - The first term is 




sinh(/3§) 


1 


dir'^Pn 8 tt^ — — 

^ n,n' X)A,A 


,7-00 27r J cosh(/?§)-I-cosh(/?(w - ^ - §)) E^(fcz)U^,(A:^-I-g^) 


X yEnikz)E^,{kz + qz) + \\’vp{kz - xb){kz + qz - X^)] + svexI^'E^{ kz){kz + qz - xb) + X{kz - xb)E^,{kz + qz)] 

„2 ,2 ^ 9 72 „2 /2 . „2 72 . „2 /2 

-=4t^ / r n—n' / \ t n —n i \ i r n—n' / B \ rn' — n / ^l.^B 


xSy -n- X'Elikz + qz))5{u: - XEl{kz))e-^{Ll7- (^)i(( -"(^) + (^)TM"( 


The second term is 

5^(7) 1 


h = 


EE 


dkz 


dw- 


sinh(/3f) 


47 r21| 7r2 27r J cosh{l3^) + cosh{l3{uj - g-^)) E^{kz)E^{kz + qz) 


2 2 / 2 7 2 \ 

xS{u; - O - X'E^kz + qz))5{u: - A£;x(fcz))AA'n'e-^Lll" (^7^ j L-"' 




After shifting the kz in the integral by bx, doing the w integral and summing over chirality, we find, 
^ sr^ sr^ f°° dkz sinh(/3§) 1 


=2 


EE 


8tt^ 2f cosh(/3f) -h cosh(/3(A£;„(fcz) - n-^)) Eykz)En'{kz + qz) 

x6{XEzi{kz) — n — X'En'{kz + qz)) ^En{kz)Eni{kz -f qz) + XX'vp{kz){kz + qz^ e 2 

. (5M) (M) (5M) ), 


2 >>■ 
(C 16 ) 


(C17) 


(CIS) 


72 = 2|(i-4^EE 


dkz 


sinh(/3§) 


47r2;| TT^ P ^,iy,J-oG 27r cosh(/3§)-|-cosh(/3(AU„(fcz) -/i - §)) E„(fcz)£’„(fcz-|-g^) 

22 ^272A ^272 

x5{XEn{kz) -n- X'En'ikz + gz))AA'n'e-^Tl-i" ( ^ U ”-”'' ^ 


(C19) 


and 
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We observe that due to the hyperbolic functions, the cooling power will be suppressed exponentially in terms of an. 
We thus focus on the contribution of the lowest Landau level to the cooling power. The only term to contribute from 
the lowest Landau level is Ii due to the vanishing of the Laguerre polynomials (L„(a;) with n < 0 is defined to be 
zero [68]). For n = 0, the first term becomes, 


L =2 


1 


■2E 


dfe. 


sinh(/3§) 


-d{-n- XvFqz)e - 2 -. 


(C20) 


Stt^ ^ 27r cosh(/3Y) + cosh{(3{XvFk^ — M ~ §)) 

After evaluating the kz integral, we arrive at our final expression for the imaginary part of the phonon self-energy, 


Im[nPh(g,L!)] = ^ 


A=± 


Stt^vfIu 


6{—rt — XvFqz)e. 2 . 


Plugging Eq. (C21) into Eq. (C26) (at zero lattice temperature), we have 


Pb = - E E / 

s,A—i q 




■5{—uj — XvFqz)e ^2 sS{uj — scsq), 


Using one of the delta functions to evaluate the w integral, we find 


Pb = - E ^ 

s,A=± 


dq f 

(2^)^ Jo 


r2TT 


d0sin0 / d(/)c^g'^A^®(scg)— 

Jo ^'Stt^vfII 


^2 02 )2 


6{—sCsq — XvFq cos 9)e 


Evaluating the (j) integral, we have 


(C21) 


(C22) 


(C23) 


pb=- y: V 

s,A=i 


dq 

(27r)2 J_, 


r-l 

„2„3 Are 


dxc‘^q'^N‘^{sCsq)- 


g^{q) .-/ . , Cs, 


-S{x + As—)e 
jf>-b 'CF 


s. 


(C24) 


Here we have taken the cutoff. A, to be the Debye wavelength, Using the last delta function to evaluate the x 
integral, we find 


P„=_ V f 


1 


dy ._ 

(2^)2^ e^^y-l 


-e 


(C25) 


^- = -E 


cVD^ 


(^fE 

r 


87r5up|4p/35(27r)2^ c 


j2 

VF \k 

) / 


2^ 


k\ 


dzz 


k+i 


- 1 


1 . 


(C26) 


Here we have used the fact that N^{—z) = —1 — N^{z) and assumed ^ -C 1. We can approximate the first term in 
the integral as 


lo e- - 1 


J dzz'‘+'‘= 2r]{k + 5)r{k + 5), 


as for large z, the integrand is small. We then have 

cVD^ 


Pb = -Y 


87r^u|.Z^4p/35(27r)2 ^ c 


oo ( \b vf \k 

(-)^E ^ (2,(fc + 5)r(fc + 5). 


fc! 




(fc + 5) ^ T, 


(C27) 


(C28) 


Assuming ^ <C 1 and performing the trivial summing over A, we have 


Pb = - 


^rr^2 rn oo ( _(b VF_T^\k 

CVD^ (!(£)5/^r5y^ 27P^ c Te) 

F'-B‘*FH \^"J C Te 


dF^Vpl‘^dpl3^{2TT)‘^ '' C ' '' Te' ^ fc!(fc + 5) 


(C29) 
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We assume that 


2 ^ 


VF Tp 


2 ^^ = 1200 Kelvin (the magnetic energy scale) 
the fact that 


7 ^ <C 1, which is a physically reasonable limit. For example, taking Tf. = 30 Kelvin, 
— = 100, and To = 140 Kelvin, we find ~ -15. Using 


<C 1, we only keep the first term in the sum. We then arrive at our expression for the power 


2^13^ C Te 

loss (after restoring factors of h) in the presence of an external magnetic field, 


Pb = - 


VD^u;% 


3207r'^u|.1^4pc^ 


(C30) 
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